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ABSTRACT 

We present a determination of the molecular gas mass in the AU Microscopii circumstellar disk. 
Direct detection of a gas component to the AU Mic disk has proven elusive, with upper limits derived 
from ultraviolet absorption line and submillimeter CO emission studies. Fluorescent emission lines of 
H2, pumped by the O VI A1032 resonance line through the C-X (1 - 1) Q{i) A1031.87 A transition, 
are detected by the Far Ultraviolet Spectroscopic Explorer. These lines are used to derive the H2 
column density associated with the AU Mic system. The derived column density is in the range 
A^(H2) = 1.9 X 10^^ - 2.8 X 10^^ cm~^, roughly two orders of magnitude lower than the upper limit 
inferred from absorption line studies. This range of column densities reflects the range of H2 excitation 
temperature consistent with the observations, r(H2) = 800 - 2000 K, derived from the presence of 
emission lines excited by O VI in the absence of those excited by Lya. Within the observational 
uncertainties, the data are consistent with the H2 gas residing in the disk. The inferred A^(H2) range 
corresponds to H2-to-dust ratios of < -^tl and a total M(H2) = 4.0 x 10"'' - 5.8 x 10~^ M^. We use 
these results to predict the intensity of the associated rovibrational emission lines of H2 at infrared 
wavelengths covered by ground-based instruments, iJS'T-NICMOS, and the Spitzer-GM&. 
Subject headings: circumstellar matter — stars: individual (HD 197481, AU Microscopii) — planetary 
systems: protoplanetary disks — ultraviolet: stars 
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1. INTRODUCTION 

Circumstellar (CS) disks around young main sequence 
stars appear to be in transition between massive, gas- 
rich protoplanetary disks and low-mass, gas-poor plan- 
etary systems. Surveys for CS disks in young stellar 
clusters suggest that gas-rich protopl anetary disks dis- 
sipate on t imescales of ~ 1 - 10 Myr (jBallv et al.lll998t 
iHaisch et a l. 2001). For solar- type stars, the timescale 
for gas dissipation is roughly equal to the theoretical time 
required for gas giant planet formation by the sta ndard 
core-accretion method (e.g. iHubickvi et al.l[2b05f ). For 
low-mass M-stars, on the other hand, the time required 
to form giant planets by core-accretion is much longer 
than it is around solar-type star s, and may be longe r 
than the typical gas disk hfetimes (jLaughhn et al.l[2004 ). 
Consequently, observations of the gas component in the 
disks of low- and intermediate-mass main sequence stars 
undergoing the transition from gas-rich protoplanetary 
disk to gas-poor debris disk are important for constrain- 
ing giant planet formation scenarios. 

AU Microscopii is a nearby {d « 10 pc) Ml star sur- 
rounded by an edge-on (inclination =1-3°) dust disk 
(jKalas et a l. 2004; Krist et al. 2005). The star is a very 
active flare star, and the majority of AU Mic studies 
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prior to imaging of it s disk focused on the stellar activity 
i)Redfield etal]l2002l . and references therein). AU Mic 
is a member of the (3 Pictoris r noving group, indicatin g 
an age of tAUMtc = 12^4 Myr (jZuckerman et al.ll2001h . 
/ 3 Pic (A5V) itself has a well-studied edge-on debris dis k 
(jRoberge et all 120001 : ILecavelier des Etangs et al.ll200H ). 
These two stars allow us to examine possible differences 
between disks of the same age around stars of very dif- 
ferent mass. 

Liu et all (|2004( l measured the 850 /im dust emis- 
sion from the AU Mic disk, and inferred a total mass 
of 0.011 Me at Tdust = 40 K. No sub-mm CO emis- 
sion was detected, implying a low CS gas mass; how- 
ever, this upper limit was not very stringent. A much 
lower upper limit on the column density of molecular 
gas was d etermined using far- UV H2 absorption spec- 
troscopy (jRoberge et al.l 120051 ). This study indicated 
that the H2-to-dust ratio in the disk is less than about 
6:1, dramatically depleted f rom t he ca, nonical interstel- 
lar ratio of 100:1. Roberge et al.l (j2005l ) mentioned that 
there was an indication of weak H2 absorption at a col- 
umn density at least an order of magnitude below their 
upper limit. They also noted the presence of very weak 
fluorescent H2 emission lines in the far-UV spectra, which 
had been seen in previous work (Redfleld et al. 2002). 

An M-star has insufficient continuum flux to ex- 
cite ("pump") far-UV fluorescent emission. In gen- 
eral, only stars of type ~B3 and earlier have the nec- 
essary spectral energy distribution (A < 1110 A) to pro- 
duce d etectabl e levels of continuum pumped UV fluores- 
cence (jFrancd [2005). However, pre-main-sequence and 
low-mass dwarf stars show high-temperature emission 
lines arising in their chromospheric and coronal regions 
llLinskv fc WxxlfT99llWood et al.l [19971 : iRedfieldeTan 
I2OO2, 2003!). These stellar emission lines can coincide 
in wavelength with transitions of H2 , providing the nec- 
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essary flux to excite detectable fluorescence when CS 
material is present. Line pumped H2 fluorescence has 
been observed from a variety of environments, including 
T Tauri disks (IValenti et al.| [2000-. Wilkinson ct al. 2003 
iHerczeg et al.ll2006l ). s unspots ([Jo rdan ct al.lll977l ). Solar 
System comet comae (iFeldman e t al. 2002|), and recently 
from planetary nebulae ~(|Lupu et al. 2006). 

Here we present a detailed description and analysis 
of the fluorescent H2 emission lines in Far Ultraviolet 
Spectroscopic Explorer (FUSE) spectra of AU Mic. The 
majority of line pumped fluorescence stud ies have fo- 
cused on H2 excitation by Lya emission 
iBlack fc van Dishoeckl Il987f ). However, in the case 
of AU Mic, we show that the observed H2 emission 
lines are pumped by the O VI A1032 stellar emission 
line, which is an important excitation source in some 
situations (IWilkinson et all l2002t iRedfield et all l2002t 
IHerczeg et al.ll2006[ ). In §2 of this paper, we briefly de- 
scribe the observations and the characteristics of the de- 
tected fluorescent H2 emission lines. In §3, we calculate 
the total fluorescent H2 flux, constrain the H2 tempera- 
ture, and determine the total column density of absorb- 
ing H2. Our estimation of the total mass of H2 gas in 
the AU Mic system appears in §4. This section also con- 
tains a discussion of the H2 heating and its implications 
for the gas location. A comparison between the observed 
AU Mic and p Pic disk properties is given in §5, and pre- 
dictions for the near- and mid-IR H2 emission line fluxes 
of AU Mic appear in §6. Our results are summarized in 
§7. 

2. FUSE OBSERVATIONS AND H2 LINE IDENTIFICATION 
FU SE performs medium-resolution (Av « 15 km s~^) 
spectroscopy in the far-UV bandpass (905 - 1187 A). The 
FUSE observatory is described in Moos et al. (2000) and 
on-orbit performance characteristics are given in Sah- 
now et al. (2000). AU Mic was observed with FUSE 
early in the mission as part of the "cool stars" pro- 
grams P118 (2000 August 26, exposure time = 17.3 
ks) and P218 (2001 October 10, exposu re time = 26.5 
ks). T hese observations are described in IRedfield et al.l 
()2002[ ). These data were acquired in TTAG mode 
through the (30" x 30") LWRS aperture. The data 
have been reprocessed with the CalFUSE pipeline v2.4.0. 
Data take n during periods of stellar flare activity were 
excluded (|Roberge et al.l [20051 ). The fluorescent emis- 
sion fines studied here fafi on the LiF 2A (1087 - 1179 
A) and LiF IB (1092 - 1187 A) detector segments. 
Archival HST Space Telescope Imaging Spectrograph 
(STIS) spectra of the chromospheric C III A1176 multi- 
plet were used to establish the wav elength calibra tion of 
the FUSE LiF IB channel (Ro berge et al.ll2005D . Due 
to a lower background in the (1 - 3) Q(3) A1119.08 A 
region, we present the LiF IB data here. 

The H2 emission lines observed in the LiF IB channel 
are shown in Figure[Tl The C~X (1 ~ 3) Q(3) A1119.077 
A and C-X (1 - 4) Q(3) A1163.807 A transitions are 
shown in the left and right panels, respectively. These 
fines are excited by absorption in the C-X (1 - 1) Q(3) 
A1031.87 A line which is coincident with the strong stellar 
O VI A1032 emission line. A description of the electronic 
ex citation of H2 and the m olecular notation are given 
in IShuU fc Beck^ithl (198^ . The (1 - 3) and (1 - 4) 
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Fig. 1.— FUSE LiF IB data for the C - X (1 - 3) Q(3) 1119.08 
A and (1 - 4) Q(3) 1163.81 A emission lines. The broad feature 
in the left panel is coronal Fe XIX ( Rcdficld et al. 2002, 2003"). 
For display, measurement uncertainties on the flux calibration are 
shown as the red error bars on the x-symbols. 



lines are expected to be the brightest ones in the O VI 
pumped fluorescent cascade (jAbgrall et al.l [l993bD . A 
Gaussian least-squares fitting routine was used to derive 
the line strength, width, velocity, and detection levels. 
The lines were marginally resolved and displayed a slight 
blueshift with respect to the stellar velocity (v* = -4.98 
± 0.02 km s"^). The (1 - 3) line was detected at 2.6 a 
with an integrated line strength of /(1-3) = 1.61 ± 0.62 
X 10-15 ergs s-i cm'^ (9.07 ± 3.49 x IQ-^ photons 
cm"^) in the (30" x 30") aperture. Its velocity'' was 
V(i_3) = -12.0 ± 2.8 km s'^ with a FWHM(i_3) = 26 ± 7 

km s~i and a peak observed fiux of 12.5 x lO^^*^ ergs 
s~ cm' -2 A-i. The (1 - 4) line was detected at 1.9 cr 
with an integrated line strength of /(i_4) = 1.48 ± 0.76 
X lO^i'^ gj.gs g-i (8.67 ± 4.45 x 10"^ photons s~'^ 

cm~2) over the same area. The (1 - 4) emission line 



was at a velocity of V(i 



(1-4) 



-22.3 ± 4.1 km s'^ with 



FWHM(i_4) = 35 ± 9 km s ^ and a peak observed fiux 

of 9.5 x 10~i5 ergs s"^ cm"^ A"^. These findings are 
summarized in Table 1. 

The observed ratios of the detected emission lines 
{R34) are consistent with the theoretical branching ratios 
within the 1 a uncertainty {R34 observed = 1.05±0.41 vs 
i?34 predicted = 0.76). The difference between the ob- 
served and theoretical ratio may be due to suppression of 
the (1 - 4) 1163.81 A fine by "the worm" seen in FUSE 
spectra from the LiF IB channel®. The line ratios de- 
rived from the LiF 2A channel are somewhat closer to 
the theoretical value, although with larger uncertainties 
(i?34 observed in LiF 2A = 0.89±0.55). The discrepancy 
is not significant within the observational uncertainties 
and does not affect the results presented in Sections 3 and 
4. The only other H2 emission line with a comparable 
branching ratio in the FUSE bandpass is the C-X (1 - 
0) (5(3) A989.73 A transition. This wavelength region 
has a lower instrumental effective area and is dominated 
by a combination of geocoronal O I and stellar N III lines. 

^ A detailed description of the FUSE 

wavele ngth calibration can be found at: 

,http : //fuse .pha. jhu. edu/analysls/calf us e .html [ 

^ More information about "the worm" can 

be found on the FUSE data analysis page: 
|http : //fuse .pha. jhu. edu/analysls/calf use_«p6 .htmll 
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TABLE 1 

H2 Emission Line Parameters from the FUSE 
Observations. 



Line Parameter 


H2 (1 - 3) Q{3) 


H2 (1 - 4) Q(3) 


Ao 


1119.077 A 


1163.807 A 


^obs 


1119.03 A 


1163.72 A 


V 


-12.0 ± 2.8 km s"! 


-22.3 ± 4.1 km 


FWHM 


26 ± 7 km 


35 ± 9 km s-l 


r a 

ilk 


9.07 ± 3.49 X 10"^ 


8.67 ± 4.45 X 10"^ 


^ peak 


7.0 X 10""' 


5.6 X 10^"' 



^ Integrated Line Strength (photons cm s Peak Line 
Strength (photons cm~^ ^~^) 



No H2 emission from C-X (1 ~ 0) (5(3) was observed. 

The H2 C-X (1 - 5) 0(3) A1208.93 A and C-X (1 - 6) 
Q(3) A1254.il A emission lines located in the STIS band- 
pass have branching ratios relative to (1 - 3) 1119.08 
A (i?53 and i?63) of roughly 0.6 and 0.1, respectively. 
Emission lines were detected at these wavelengths at the 
2 ~ 3 cr level. Their velocities and line strengths are con- 
sistent with the observed H2 emission lines in the FUSE 
spectra. We will focus on the lines in the FUSE band- 
pass for the remainder of the paper, but we note that 
the detection of additional H2 fluorescent emission lines 
with an independent instrument makes the conclusions 
presented in Sections 3 and 4 more robust. 

3. ANALYSIS AND RESULTS 

3.1. Total Emitted H2 Flux 

Using the measured H2 emission line fluxes, we can cal- 
culate the total fluorescent output from the O VI pumped 
cascade. The total emitted flux out of the electrovibra- 
tional state {n',v',J'), J2 ^ij^ is given by 
3 

where i refers to the upper state (n',v',J'). The indices 
j, k, and I refer to the lower states {n" ,v" ,J"). lik is 
the measured flux of the FUSE band lines (in photons 
s""'^ cm~^ ). The ratios of ind ividual to total Einstein 
A-values (jAbgrall et al.lll993bf) are the branching ratios, 
and is a correction for the e fficiency of predissocia tion 
in the excited electronic state ()Liu fc DalgarnolllQQGf ) . In 
the case of AU Mic, we are concerned with Werner band 
emission {n' - n" ^ C - X = C^U^ - X^Y.+), v' ^ 1, 

and v" = 3 and 4 for the 1119.08 and 1163.81 A lines, 
respectively. The predissociation fraction for the Werner 
bands is zero (fc — 0; Ajello et al., 1984). 

Following this procedure, we arrived at the total emit- 
ted photon flux, ^ /ij, derived from the observed (1 - 3) 
3 

1119.08 A and (1 - 4) 1163.81 A emission lines individu- 
ally (4.3 ± 1.7 X 10-4 and 3.1 ± 1.6 x lO""* photons s'^ 
cm for the (1 - 3) and (1 - 4) lines, respectively.) For 
the determination of the total H2 column density pre- 
sented below, we take the average of these values, ^ lij 

3 

— 3.7 ± 2.3 X 10^"' photons s~^ cm^^. The error on 
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Fig. 2. — Excitation of fluorescent H2 emission lines by stellar 
O VI emission. The intrinsic stellar O VI A1032 emission model 
(Roberge et al., 2005) is shown by the solid blue line. The H2 ab- 
sorption model with a column density in the {v,J) = (1,3) level of 
Nx,i,3 = 3.48l^;4| X 10^3 cm-2 is shown at the top with the black 
dot-dash line; this model has been convolved with the FUSE LSF. 
The solid red line shows the stellar emission model times the H2 
absorption model; the shape of the stellar emission line is modified 
by superimposed absorption from the H2 C - X {1 — 1) Q{3) tran- 
sition at 1031.87 A. This absorption line is the excitation channel 
for the observed fiuorescent H2 emission lines shown in Figure 1. 
The red dashed lines show total models with absorption columns 
corresponding to the quoted error ranges on the H2 column density 

the total emitted flux was determined such that both 
the (1 - 3) 1119.08 A and (1 - 4) 1163.81 A values are 
consistent with the mean value; this approach ensures a 
conservative estimate of the measurement uncertainties. 

3.2. Inferred H2 Column Density 

In order to determine the total column density that is 
associated with the observed level of emission, we made 
three assumptions: 1) that fluorescence is the only source 
of the observed emission, 2) that the stellar O VI A1032 
emission line is the only source of pumping photons, and 
3) that the total number of photons is conserved. The 
first assumption (1) seems warranted as H2 cannot be 
electronically excited by shocks or collisional processes 
with other gas or dust particles. Additionally, excitation 
by elect ron collisions has a distinct far-UV emission sig- 
nature (lAjello et al.lll982( ) which is not observed towards 
AU Mic. (2) AU Mic does not emit stellar continuum 
at wavelengths coincident with the absorbing transitions 
that produce the fluorescent emission lines, supporting 
the assumption that the O VI emission line is responsi- 
ble for the observed excitation. (3) Detailed calculations 
of the formation and destruction of molecules in the AU 
Mic system are beyond the scope of this work, thus we 
a ssumed photon conser vation for what follows. 

iRoberge et all (|2005[ ) modeled the O VI AA1032/1038 
and C II AA1036/1037 stellar emission lines in the FUSE 
data by fitting a combination of narrow and broad Gaus- 
sians to each line. The high-excitation ionic lines thought 
to originate in the chromosphere and transition regions 
of low mass dwarf stars are known to be poorly fit by 
a single Gaussian component (Linsky & Wood, 1994; 
Redfield et al., 2002). The narrow component of the 
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profile has a FWUMnarrow = 44 ± 6 km s^^ with 
a velocity of v„arrow — -4.6 ± 1.6 km s^^. This is 
consistent with the stellar velocity (v, = -4.98 ± 0.02 
km s~^). The broad component has FWUMbroad — 
109 ± 25 km with a some what redshifted veloc - 
ity, Vbroad = +7.7 ± 9.3 km s'^ (jRoberge et al.ll2005[ ). 
This two-component model was used as the stellar emis- 
sion profile to be absorbed by the C-X (1 - 1) Q{3) 
A1031.87 A transition. The absorption pro file was cre- 
ated from an H200IS optical depth template (jMcCandlissI 
12003) with a conservative Doppler pa rameter 6 = 2 km 



s ^ (jLecavelier des Etangs et al.ll2001[ ). TheH20olstem 



plates are optical depth arrays that can be used to fit ar- 
bitrary H2 column densities for many rovibrational states 
for 6- values from 2 - 20 km s'^ (|McCandhssll2003f ) . The 
stellar emission model was binned to 0.01 A pixels to 
match the H200IS wavelength grid. A grid of column 
densities was searched to find the minimum difference 
in the absorbed [1(0 Ylmodei) ^ 1(0 Vlmodei x H2 ab- 
sorption)] and emitted fluxes (the equilibrium condition). 
This method found N{R2,n,v,J) = Nx,i,3 = i-^^tlH 
X 10^'^ cm^^. The O VI model and H2 absorption are 
illustrated in Figure [2] 

The total A^(H2) can be determined from Nx,i,3 by 
assuming a Boltzm ann distribution and an estimate for 
the temperature. iLiu et ahl ([2004) assumed that the 
CO in the AU Mic d isk is in the rmal equilibrium with 
the dust at 40 K and lRoberge et al.l (|2005l) consider ab- 
sorption out of the i; = 0, J = 0, 1 states of H2 assuming 
that r(H2) < 200K. To make our own determination of 
the H2 excitatio n temperature, we u sed far-UV H2 fluo- 
rescence models (iFrance et al.ll2005l ) to predict the tem- 
perature dependence of the H2 emission spectrum. The 
relative fluxes of the emission lines are determined by the 
shape and strength of the exciting radiation field and the 
H2 abundances in the rovibrational levels of the ground 
electronic state. The level populations are determined 
by the H2 column density and excitation temperature. 
We computed fluorescence models for a range of excita- 
tion temperatures (40 < T(H2) < 2000 K) and column 
densities (10^^ < A^(H2) < 2 x 10^^). The fluorescence 
code used the 1030 - 1040 A O Vl+C II stellar emission 
model, described above, as the exciting radiation field. 

At excitation temperatures < 700 K, fiuorescent emis- 
sion lines excited by absorption out of the v — and J 
= 0, 1, and 2 levels dominate the output in the 1100 - 
1187 A wavelength range. These lines are not seen in the 
FUSE spectra. The observed (1 - 3) and (1 - 4) emis- 
sion lines become the strongest at temperatures above 
800 K, providing a lower limit on the excitation tempera- 
ture. This spectral variation with excitation temperature 
is due to the distribution of higher rovibrational levels 
within the ground electronic state. Once the excitation 
temperature is high enough to significantly populate the 
(u,J) — (1,3) level, this fiuorescent route dominates due 
to the coincidence with O VI A1032. An upper limit on 
T(H2) can also be set from the observed spectral charac- 
teristics. The lack of Lya pumped fluorescence indicates 
(§4.2) that the region producing H2 line emission is cooler 
than « 2000 K. 

In order to present a fiducial column density value 
that gives a sense of the measurement errors, we adopted 
T(H2) — 1000 K as a characteristic temperature, which 



gives a total 7V(H2) =4.24+^:™ x lO^^ cm^^. We empha- 
size however, [T(H2),A^(H2)] combinations ranging from 
[800 K, 1.93 X 10" cm-2] to [2000 K, 2.80 x lO^^ cm^^] 
are consistent with the data (see §4.2 for a more detailed 
discussion of the molecular gas temperature). Column 
density distributions are shown in Figure 3 for the 800 
and 2000 K cases. 

4. PHYSICAL PROPERTIES OF THE MOLECULAR GAS 
COMPONENT 

4.1. Mass 

A dust mass (Md) of 0.011 M m in the AU M ic disk 
(at d < 70 AU) was measured bv lLiu et al.l (|2004l ) using 
850 ^m SCUBA observations. More recently, the dust 
mass in the AU Mic disk has been estimated from the 
visible and near-IR scattered light profiles of the disk. 
Calculations based on scattered light find masses smaller 
(~ a few - 70 X 10"'* depending on the grain prop- 
erties and size distribution; Augereau & Beust 2006) or 
equal (~ 0.01 M®; Strubbe & Chiang 2006) to the sub- 
mm value. In o rder to dire ctly compare with H2 ab- 
sorption studies (jRoberge et al. 200 5), we will u se 0.011 
as the CS dust mass here. Liu et al.l ()2004f ) also set 
an upper limit on the CO column density in the disk of 
A^(CO) < 6.3 X 10" cm-2 from a CO (3 ~ 2) 346 GHz 
emission non-detection. Assuming an iV(C0)/iV(H2) ra- 
tio of 10~^, they place an upper limit on the H2 col- 
umn in the disk of A^(Il2) < 6.3 x 10^° cm-^. Their 
assumption is supported by recent studies of the dif- 
fuse ISM at comparable values of A^(CO) where the 
A^(C0)/A^(H 9) ratio is ob served to be in the range of a 
few X 10-^ (Burgh et al.ll2007f) . iLiu et all (|2004[ ) place 
an upper limit on the mass of H2 gas in the disk of 
< 1.3 Mq and the H2-to-dust ratio in the disk at 
M/z^/Mci < 118:1. The limit on this ratio was further 
decreased by II2 absorption line spectroscopy (Roberge 
et al. 2005; A^(H2) < 1.7 x lO^^ cm'^) to MnJMd < 
6:1. 

Comparing the result derived in §3.2 with that 
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Fig. 3.— The H2 column density distributions for T(H2) = 800 
(diamonds) and 2000 K (squares). Excitation temperatures be- 
low 800 K are inconsistent with the observed fiuorescent spectrum. 
Temperatures greater than 2000 K are ruled out by the absence 
of H2 emission lines pumped by Lya. Error bars on the column 
density are shown for the {v,J) = (1,3) level. 
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TABLE 2 

H2-T0-DUST Ratio Determinations for the AU Mic System. 



M(H2) 


Technique 




H2-to-Dust 


Reference 






(Me) 






< 1.3 


CO (3 - 2) Emission 


0.011 


< 118:1 


Liu et al., 2004 


< 0.07 


H2 UV Absorption 


0.011 


< 6:1 


Roberge et al., 2005 


4.0 X 10"* 


H2 Fluorescence 


0.011 


0.036:1 


This Work, T(H2) = 800 K 


5.8 X 10"® 


H2 Fluorescence 


0.011 


5.2 X 10-'*:1 


This Work, T(H2) = 2000 K 



from 850 fj.ni dust emission (|Liu et al.ll200i ) 



of iLiu et al.l (|2004D . we infer values for 7V(H2) that 
are « 3.3 x 10^ - 2.3 x 10^ below their upper hmit. 
The corresponding H2-to-dust ratios are Mj^^/A^d = 
(0.036 - 5.2 X lO"*):!, or MnJMd < -^-.l. This 
gives a total mass range of M(H2) = 4.0 x 10"''' - 
5.8 X 10-6 Me- The value for r(H2) = 1000 K is M(H2) 
= S.T^'^ X 10^^ M0. These results are summarized in 
Tabled' 

4.2. Temperature 

Observations of H2 emission excited by Lya through 
the 5 - X (1 - 2) P(5) 1216.07 A and B - X (1 - 2) R{6) 
1215.73 A coincidences are generally thought to indicate 
an H2 ground state populatio n characterized by tem- 
peratures rfHz ) > 2000 K (Black fc van Dishoecklll987l: 
Wood fc Karovs ka 2004i: iHerczeg et al.ll2004l : iLupu et all 



20061) . In the case where O VI and Lya excitation are 



both observed, Lya excitation usually dominates (e.g.- 
T Tauri stars; Wilkinson et al. 2002; Herczeg et al. 
2005). The observation of O VI pumped fluorescence in 
conjunction with the absence of Lya fluorescence allows 
us to constrain the molecular gas temperature in the AU 
Mic disk. As discussed in §3.2, the observed fluorescence 
spectrum sets the lower limit on T(H2) > 800 K, and we 
suggest that the lack of Lya pumped lines in the FUSE 
data imply an upper limit of r(H2) < 2000 K. 

A quantitative calculation of the expected flux from 
Lya induced fluorescence is complicated b y strong inter- 
stellar H I absorption of the line profile (jPagano et al.l 
I2OOOI : iRedfi cld ct al. 2002,). It seems clear that the local 
AU Mic Lya radiation field is at least an order of mag- 
nitude more intense than the local O VI radiation field 
(Figure 5 in both Pagano et al. 2000 and Redfield et 
al. 2002). The A- values for the strongest Lya pumped 
lines a re similar (to within 30 %) to those pumped by 
O VI (lAbgrafl et all Il993a| bl). Ignoring optical depth 
and extinction effects, the column densities in the absorb- 
ing transitions control the resultant emission spectrum. 
We can define the ratio of column densities in the rele- 
N{1,3)/{N{2,5) + N{2,6)). When 

< 10 (this value is set by the rough estimate of 
stellar Lya-to-0 VI), we would expect a detectable con- 
tribution from Lya pumped fiuorescence. For gas tem- 
peratures of T(H2) = [1000, 2000, 3000, 4000 K], the 
corresponding ratios are R^yJ = [545.4, 16.0, 4.9, 2.7]. 
Thus, only when T(Il2) < 2000 K do we expect to detect 
II2 emission from O VI excitation in the absence of lines 
excited by Lya. 

4.3. Spatial Origin and Heating 



vant states R^^^ 

^Lya 



OVI Model - 
OVI Model • H, Template - 
1-0 Errors on N(H,) 



xFUSE Data 




1031.5 1031.8 1032.0 _ 1032.2 1032.4 

Wavelength (A) 

Fig. 4. — A comparison between the predicted and observed 
O VI line profiles. The blue and red lines are the same as in 
Figure 2. The FUSE data are overplotted with x's. The required 
H2 absorption is consistent with the observed line profile, within 
the measurement errors. This indicates that the H2 absorption 
could arise completely within the AU Mic CS disk. 



The H2 required to produce the observed fluorescence 
need not be coincident with the stellar line of sight. The 
emitting gas could reside in a cloud that extends beyond 
the disk. If the necessary column density suggests ab- 
sorption that is not observed along the line of sight, this 
could be evidence for an extra-planar gas component. In 
order to explore this possibility, we compared a model 
of the O VI line proflle modified by the required H2 ab- 
sorption to the O VI profile observed in the FUSE spectra 
(shown in Figure 4). This assumes that the total required 
absorbing column density lies along the line of sight, as 
it would if the II2 gas is entirely in the disk. The O VI 
model with the required H2 absorption superimposed is 
consistent with the data, within the measurement uncer- 
tainties. This leaves open the possibility that all of the 
emitting H2 gas lies in the disk, and only the relatively 
small absorbing column density and low signal-to- noise of 
the FUSE spectra prevent it from being clearly detected 
in absorption against the O VI 103 2 emission line, as ha s 
been observed in other CS disks (jRoberge et al.l[200l[ l. 
It is interesting to note that the column density range 
we derive here (7V(H2) = 1.9 x lO" - 2.8 x lO^^ cm^^) 
is consistent with the pos sible line of sight ab sorption 
suggested in Section 4.2 of lRoberge et al.l (|2005l ). 

We can also test if our T(H2) estimate is consistent 
with a disk origin for the emitting gas. r(H2) ^ 1000 
K is warmer than typical debris disk gas temperatures 
considered in previous theoretical {T < 300K, Kamp 
& van Zadelhoff 2001) and observational (T - 40 K, 
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Liu et al. 2004; < 200 K, Roberge et al. 2005) stud- 
ies. The first possibility is that the gas is in close prox- 
imity to the star. Calculations have been presented 
for the detectability of Lyg pumped H2 i n the spec- 
tra of late- type stars bv iJordan et al.l (|1978l ). but as we 
have shown, the temperature of the molecular gas as- 
sociated with AU Mic is below the level required for 
that process. Starspots have lower effective tempera- 
tures than the conventional photosphere, and the AU Mic 
spot temperature has be en measured in variability stud- 
ies ()Rodono et al.l [l986l ). However, the observed spot 
temperature, Tgpot = 2650 K, is still substantially above 
our 2000 K limit. The ratio of spot-to- unspotted tem- 
peratures in AU Mic is ~ 0.76. This value agrees well 
with the spot-to-unspotted temperature ratios found for 
more massive active stars (~ 0.66 ~ 0.86; NefF et al. 1995; 
O'Neal et al. 1996) from an analysis of TiO absorption 
bands. We find that even in the coolest regions of the 
AU Mic surface, a photospheric origin for the observed 
H2 emission can be most likely ruled out. 

The AU Mic corona is characterized by electron tem- 
perat ures in the range of lO'' < < 10^ K (jMaran et al.l 
Il994[ ). Electron densities in the transition region 
and c oronal regime ar e of order rie = 5 x 10^" 
cm~^ (jMaran et al.lll99"^ . with suggestions of values sev- 
eral orders of magnitude greater (Redfield et al. 2002). 
We can estimate the characteristic survival time for 
molecules near the AU Mic transition and coronal regions 
by assuming that the H2 is at low densities (i.e.- not con- 
tained in dense clumps), using e~ + H2 cros s sections to 
deter mine the collisional dissociation rate tMarti iTet al.l 
Il998f) . Taking a conservative value for the electron tem- 
perature (10'* K), we find the e-folding time for H2 dis- 
sociation via electron impact, Tdiss = (7e'T-e)~^, where je 
is the electron impact dissociation rate coefficient, to be 
Tdiss ^ 1 minute. It is possible that the interstellar condi- 
tions assumed to estimate this time scale do not apply di- 
rectly to the AU Mic environment, however the basic pic- 
ture that molecules cannot survive in the immediate envi- 
ronment of AU Mic seems robust even if the parameters 
vary by several orders of magnitude. As an additional 
constraint, we consider the possible observational conse- 
quences assuming that molecules could survive in regions 
as close as a few stellar radii. H2 heated by collisions 
in the stellar atmosphere might be expected to show the 
well charac terized signature of electron-irtipact i nduced 
excitation (jAiello et all 119821 : lAjello et al.l ll98l ). The 
electron-impact excitation spe ctrum of H2 peaks in the 
FUSE and STIS bandpasses (|Gustin et al.lf200l l200l . 
and no emission from these features is observed. 

Some gra in species may be resistant to sublimation at 
ann et al. Il200g). and it may be possible 
that the H2 is heated by collisions with grains in thermal 
equilibrium with the stellar radiation field (r « 6 i?* for 
T — lOOOK). However, previous studies have found no 
dust emission or scattered light in the inner disk near 
AU Mic. In the optica l, the disk is cleared inside r < 7.5 
AU (|Krist et al.ll2005D . Probing the 850 ^ni dust emis- 
sion, the inferred inner radius is r < 17 AU (|Liu et al.l 
[2001 . 

The inferred range of r(H2) could also be produced 
by heating processes operating in the disk. There have 
been considerable efforts towards modeling the gas and 
dust components of CS disks (Kamp & van Zadelhoff 



2001; Jonkheid et al. 2004; Besla & Wu 2007; and 
references therein), although most of these efforts have 
focused on higher mass stars (Herbig Ae and T Tauri 
stars). Nevertheless, we will use these models to under- 
stand the important processes heating the molecular gas 
in AU Mic, noting where certain assumptions are invalid 
for an M star disk. Heating processes include photo- 
electric heating by grains, collisional de-excitation of H2, 
photodissociation of H2, H2 formation, gas-grain colli- 
sions, ca rbon io nization, gas-grain drift, and cosmic ray 
heating ()Kamp fc van Zadclho fi 2001) . Relevant cooling 
processes that regulate the gas temperature include [O I] 
cooling (from A6300 A and A63.2 iim), [C I] and [C II] 
cooling, H2 rovibrational line cooliiig, Ly g cooling, and 
CO cooling (Kamp fc van Zadelhofaf200l[ ). 

Gas-grain collisions only heat the gas when Tgas < 
Tdust, and there is uncertainty whether gas-grain d rift is 
an efficient heating mechanism ()Besla fc Wull2007l ). The 
gas temperature would have to be much higher than we 
observe for Lya or [O I] A6300 A cooling to contribute 
significantly, and the lifetimes of the rovibrational states 
of H2 are very long, making these IR transiti ons ineffi- 
cient coolants. CO emission is not detected ()Liu et al.l 
[200l . so it is hard to assess this contribution to the 
cooling. Finally, AU Mic lacks the far-UV (912 - 1110 A) 
stellar continuum that drives the photodissociation of H2 
and CO, and produces C II through the photoionization 
of carbon. Processes that depend on this flux [heating : 
photodissociation of H2 and carbon ionization, cooling : 
[C II] emission) must be driv en solely by the interstellar 
radiation field (lDrainel[l978D . Relative to models of A- 
star disks, we presume these processes are of diminished 
importance in AU Mic. 

This leaves photoelectric and H2 formation heating, 
and far-IR fine structure line cooling as the dominant 
processes that determine the gas temperature in the AU 
Mic disk. This co nclusion generally a grees with the sce- 
nario put forth bv lBesla fc Wul ()2007f ). who calculate gas 
temperatures in the disks of more massive stars (K2 and 
earlie r) of Tppg > 400 K, dep ending on the model param- 
eters. iJonkheid et al.l (|2004f) find that strong photoelec- 
tric heating leads to gas temperatures as high as 1000 
K in the surface regions of disks, although they consider 
more strongly fiared disks than AU Mic. We note that 
1000 K is similar to the temperature found for interstel- 
lar H2 where grain-formation pumping is thought to b e 
a dominant excitation source (jSpitzer fc Cochranlll973f) . 
Suffice to say, photoelectric and H2 formation heating 
seem to be capable of elevating the gas temperature to 
roughly the observed level, but more modeling work is 
needed for disks around low mass stars. Combining this 
with the arguments against a stellar origin for the ob- 
served H2 emission given above, we favor the hypothesis 
that the fiuorescent emission originates in the CS disk. 

5. COMPARISON WITH /3 PIC 

AU Mic is a member o f the (3 Pictoris moving 
group (|Zuckerman et al.ll200in . meaning the AU Mic disk 
is roughly the same age as the well-studied (3 Pic debris 
disk. (3 Pic ( A5V) is roughly 3 .6 times more massive 
than AU Mic (jKalas et al.l l2004D . In this section, we 
briefly compare the molecular gas properties derived for 
AU Mic with previous observations of gas in the (3 Pic 
CS disk. 
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TABLE 3 

Predicted Near and Mid-IR H2 AU Mic Emission Line Brightnesses. 





Line 


Wavelength 


Instrument 




Brightness^ 




{v' 


- «") S{.J) 


(/im) 




(ergf 


i cm" s" ^ sr~ 


') 


(1 


- 0) 5(1) 


2.12 


HST-mCMOS /Ground 




9.01 X 10"^ 




(2 


- 1) S(l) 


2.25 


Ground 


1.13 X 


10-9 - 7.47 X 


10-8 


(0 


- 0) 5(7) 


5.51 


Spitzer-IRS 


5.74 X 


10-7 - 4.93 X 


10-7 


(0 


- 0) 5(6) 


6.11 


Spitzer-IRS 


2.95 X 


10-7 - 1.48 X 


10-7 


(0 


- 0) 5(5) 


6.91 


Spitzer-IRS 


2.02 X 


10-'5 - 3.30 X 


10-7 


(0 


- 0) 5(4) 


8.03 


Spitzer-IRS / Ground 


8.84 X 


10-7 - 6.40 X 


10-8 


(0 


- 0) 5(3) 


9.66 


Spitzer-IRS 


2.29 X 


10-'3 - 8.13 X 


10-8 


(0 


- 0) 5(2) 


12.28 


Spitzer-IRS / Ground 


3.84 X 


10-7 - 7.39 X 


10-9 


(0 


- 0) 5(1) 


17.03 


Spitzer-IRS / Ground 


2.55 X 


10-7 - 2.99 X 


10-9 


(0 


- 0) 5(0) 


28.22 


Spitzer-IRS 


4.25 X 


10-9 - 3.41 X 


10-11 



^ Predicted brightness ranges reflect the adopted temperature 
X IQi^ - 2.8 X IQi^ cm"^ described in §3.2. 

Multi-wavelength model fits to the Pic CS disk SED 
predict a dust mass of 0.037 (jPent et al.ll2000l ). The 
molecular hydrogen mass in P Pic is less clearly de- 
fined. ISO observations of the mid-IR emission lines 
of H2 find a large molecular gas reservoir (M(H 2) ~ 
57 Me) associated with /3 Pic (|Thi et all l200i[ ). It 
seems unlikely that this emission is distributed uniformly 
throughout the disk. Using UV absorpt ion techniques 
analo g ous to those presented for AU Mic (jRoberge et al.l 
I2OO5D , iLecavelier des Etangs et all (|200lD report a non- 
detection of H2 absorption in the edge-on disk. They set 
an upper limit on the H2 column density of A^(H2) < 10^* 
cm~^, corresponding to a molecular gas mass of < 0.095 
Me. 

The ratio of the /3 Pic and AU Mic dust masses (~ 3.4) 
is approximately equal to the ratio of their stellar masses 
3.6). The observed dust in both disks is thought to be 
continually replenished. Collisions of larger parent bod- 
ies in the disk can rcpopulate the small grain population 
that is detected as scattered light in the optical/near-IR 
and as thermal emission at longer wavelengths. Models 
have shown that this scenario can reproduce the observed 
dust prop erties of Pic ( Thcbault et al. 200 ^ and 
AU Mic (jAugereau fc BeustI 120061 : IStrubbe fc Chiangj 
|2006[) . Grain collisions may also be responsible for 
reple nishing the metallic gas observed in the /? Pic 
disk (jFernandez et al.ll2006f) . but it is unclear if this pro- 
cess contributes to the gas phase H2 abundance. 

The H2-to-dust ratio in the p Pic disk (from the UV 
absorption line upper limit) is < 3:1. This is consistent 
with our range of H2-to-dust ratios for AU Mic (0.036 - 
5.2 X IQ-^:!). If H2 gas is present in the /3 Pic disk at 
a similar gas-to-dust ratio as AU Mic, a natural ques- 
tion arises: Why were the fluorescent emission lines not 
detected in FU SE observations of (3 Pic? The answer 
is most likely related to the earlier spectral type of (3 
Pic. The stellar O VI emission line from an active M- 
star such as AU Mic is considerably stronger than in /3 
Pic. The peak flux at the O VI A1032 A line center (co- 
incident with the absorbing II2 transition studied here) 
is over an order of magnitude higher in AU Mic. Ad- 
ditionally, t he stellar continuum of Pic ext ends down 
to 1100 A (jLecavelier des Etangs et aT1l2001[ ). lowering 
the line-to-continuum ratio at the strongest emission line 
wavelengths. 



iges of r(H2) = 800 - 2000 K, corresponding to iV(H2) = 1.9 

6. IR BRIGHTNESS PREDICTIONS 

H2 does not have an intrinsic dipole moment, hence the 
rovibrational transitions of the molecule proceed by the 
slower quadrupole channel, making t hem optically thin in 
most astrophysical environments fjBlack fc van Dishoeckl 
Il987f) . Assuming the optically thin case, we used the 
derived column density distributions ([T(H2),A^(H2)] = 
[800 K, 1.93 X 10" cm-2] - [2000 K, 2.80 x 10^^ 
cm^^]) to predict the near and mid-IR II2 emission 
line strengths for transitio ns that a rc observabl e from 
ground-based facilities ( Speck et al.l 2003: AUer s et al.l 
I2005D , HST-mCMOS (iMeixner et al.ll2005D . or Spitzer- 
IRS (iHouck et all 120041: iHora et all \20o4 . The most 
readily observable H2 lines are the rovibrational lines (1 - 
0) S{1) A2.12 ^m and (2 - 1) S{1) A2.25 ^J.m, and the 
pure rotational lines (0 - 0) S'(7) - S'(O) A5 - 29 fim. The 
predicted line intensity can be calculated from 

/ Npij"+2) Ag \ he 

imsij")^(^ — a; — ji^ 

where ImSjJ") is in unit s of ergs s~^ cm~^ 
sr-^ ([Black fc van DishoecS 119871 : [Rosenthal et all 
[2OOOI) . m is the {v' - v") transition, p refers to the 
upper vibrational level, and q labels the A-value and 
wavelength for the relevant transition. The J" + 2 
notation is a consequence of the S br anch (A J = -1-2) 
transi tion. The A-values are from IWolniewicz et all 
(|1998f ). The predicted line strengths are given in Table 
3. 

The exact detection limits of these lines will be de- 
termined by the angular size of the AU Mic emission, a 
larger filling fraction will increase the observed signal at 
a given surface brightness. The Spitzer-lUS has access 
to the largest number of these lines. Even assuming op- 
timistic aperture filling fractions of unity, the brightest 
of these lines are predicted to have Y^, < 0.05 mJy, which 
is at or below the IRS noise limit. This result is consis- 
tent with the non-detection of H2 emission from the 15 
CS disks around young Sun-like stars in the Formation 
and E yolution of Planetary Systems Spitzer Legacy Pro- 
gram l|Pascucci et al.ll2006[ ). iJS'T-NICMOS has the ca- 
pability for narrow band imaging in the (1 - 0) S'(l) A2.12 
/im line, however we predict that achieving contrast with 
the stellar emission in this bright star [mv = 8.8; Kalas 
et al. 2004) will be difficult. At a temperature of Te// = 
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3500 K, the photosphere wiU emit strongly in the near- 
IR (X„iax = 0.83 /im). Near-IR imaging will require high 
dynamic range in flux to achieve contrast between the 
stellar emission and the faint molecular gas. Calcula- 
tions for other instruments can be performed using the 
values in Table 3. We suggest that high resolution near- 
IR spectroscopy could be the most promising technique 
for future observations of H2 in the disk. 

7. SUMMARY 

We have presented far-UV observations of H2 emission 
in the AU Microscopii CS disk. The spectra displayed 
fluorescent emission lines excited by stellar O VI A 1032 
photons coincident with the C-X (1 - 1) Q(3) A1031.87 
A transition. These lines imply a total column density 
in the molecular gas of N(H2) = 1.9 x lO" - 2.8 x 10^^ 
cm~^. This detection is roughly two orders of magnitude 
smaller than published upper limits on N(H.2). Compar- 
ing this value with previous limits on the gas mass in the 
system, we find M(H2) = 4.0 x 10"^ - 5.8 x 10"^ M®. 
Using the molecular mass and the 850 /xm dust emis- 
sion, we found a gas-to-dust ratio of < ^:1. The derived 
column densities and gas masses depend upon the as- 
sumed excitation temperature, which we estimate to be 
800 - 2000 K. We presented the basis for this tempera- 



ture distribution, and discussed the value in the context 
of the AU Mic system. We conclude that the warm H2 is 
most likely associated with the disk, with photoelectric 
heating and formation pumping as the dominant heat- 
ing mechanisms. The intensity was predicted for several 
near and mid-IR lines of H2 . These lines are accessible to 
current ground and space-based observatories, although 
due to the low column density and weak intrinsic nature 
of the lines, they will be challenging to detect. 
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